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Mode-coupling as a Landau theory of the glass transition 
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PACS 64 . 70 . Q — Theory and modeling of the glass transition 
PACS 64 . 60 . Ht - Dynamic critical phenomena 

PACS 64 . 60 . F — Equilibrium properties near critical points, critical exponents 

Abstract. - We derive the Mode Coupling Theory (MCT) of the glass transition as a Landau 
theory, formulated as an expansion of the exact dynamical equations in the difference between the 
correlation function and its plateau value. This sheds light on the universality of MCT predictions. 
While our expansion generates higher order non-local corrections that modify the standard MCT 
equations, we find that the square root singularity of the order parameter, the scaling function in 
the (3 regime and the functional relation between the exponents defining the a and /3 timescales 
are universal and left intact by these corrections. 
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The Mode-Coupling Theory of glasses (MCT), devel- 
oped since the mid-eighties following the seminal work of 
Gotze [1] and Leutheusser [2], has significantly contributed 
to our understanding of the slowing down of supercooled 
liquids. One of its cardinal predictions is the appearance 
of a non trivial /3-relaxation regime where dynamical cor- 
relation functions pause around a plateau value before fi- 
nally relaxing to zero. In the vicinity of this plateau value, 
the theory predicts two power-law regimes in time (or in 
frequency), and the divergence of two distinct relaxation 
times, r a and 773, at the MCT critical temperature T^. Al- 
though this divergence is smeared out by activated events 
in real liquids, the two-step relaxation picture suggested 
by MCT seems to account quite well for experimental and 
numerical observations [3], at least in weakly supercooled 
liquids and for hard sphere colloidal systems. 

Originally, MCT was obtained as an uncontrolled self- 
consistent approximation within the Mori-Zwanzig pro- 
jection operator formalism for Newtonian particles. This 
scheme yields an integro-differential equation for the dy- 
namic structure factor C(k, t) that captures mathemati- 
cally the slowing down of the dynamics and the appear- 
ance of a two-step relaxation at equilibrium. It pro- 
vides very detailed predictions for the scaling properties 
of C(k, t) in the vicinity of the plateau value /kSk, where 
S*k is the static structure factor and /k is called the non- 
ergodicity parameter (akin to the Edwards- Anderson pa- 



rameter in spin glasses) . Alternative derivations of MCT 
based on field theory have been sought for [4] and research 
on this topic has continued until now [5-8] . It was also re- 
alized that the same integro-differential equations describe 
the exact evolution of the correlation function of mean- 
field p-spin glasses [9]. This is important for at least two 
reasons: 

• Technically, it shows that the MCT approximation is 
realizable: there is a well defined system for which it 
is exact; hence MCT does not violate basic physical 
constraints. 

• Physically, it brings in a very useful interpretation 
of the MCT freezing transition in the "energy land- 
scape" parlance. Above the transition X^, the dy- 
namics is dominated by unstable saddles that become 
progressively less and less unstable as one approaches 
the transition; below the transition, there are only lo- 
cal minima that are separated by infinite barriers in 
mean-field, so that the system is forever trapped in 
one of them [9]. For non mean-field systems, these 
barriers are finite and the transition is smoothed. 
MCT can be naturally embedded within the broader 
Random First Order Theory [10] of the glass transi- 
tion. In this context it describes the high temperature 
region where metastable states are still in embryo. 

However, this analogy shows that MCT is (at best) an 
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incomplete theory of real supercooled liquids, and needs 
to be corrected and enhanced. An important question is to 
know to what extent the quantitative predictions of MCT 
are stable with respect to the ignored contributions. In 
fact, it is genrally accepted that the quantitative value of 
the critical temperature (or the critical density) obtained 
solving MCT equations is incorrect, as is the interaction 
parameter A, which is a functional of the static structure 
factor and fixes the value of critical exponents. However, 
when MCT is used to fit empirical data, it is assumed 
(without much justification) that the predictions about 
the critical behavior remain correct if Td and A are treated 
as adjustable parameters. This implicitly assumes that 
some MCT predictions are universal, e.g. the square root 
singularity of /k and the relation between the exponents 
describing the divergence of r a and t@, and others which 
are not, e.g. Td and the actual value of the exponents! 

Clearly, the present theoretical understanding of MCT 
needs to be improved. Assessing the degree of structural 
stability of the theory and its universality properties is a 
crucial issue to resolve both for theoretical and practical 
purposes. One step in this direction has been performed 
by Szamcl [11] and then later generalized by Mayer et al., 
where a schematic version of MCT including higher order 
correlations was proposed and analyzed [12]. The result is 
that whenever the theory is truncated at any finite order in 
the n-boby correlations, the phenomenology of MCT is ex- 
actly recovered with a finite Td, whereas Td = when the 
theory is treated exactly to all orders. Another approach 
was followed in [6], based on a field-theoretical formula- 
tion of MCT consistent with the Fluctuation-Dissipation 
Theorem, which suggests closure schemes different from 
standard MCT. 

The aim of this work is to argue that MCT can be 
rephrased as a Landau theory of the glass transition, based 
on general assumptions about the nature of the dynam- 
ical arrest but without relying on any particular model. 
Therefore, some predictions of MCT are indeed generic 
and should be useful in a certain regime of time and tem- 
perature. Note that a Landau theory for the glass tran- 
sition has been developed also in [13], but it has a very 
different starting point and it does not focus on MCT. 

The Landau theory is a general phenomenological ap- 
proach to equilibrium phase transitions [14, 15]. It re- 
lies on a number of natural hypotheses, such as symme- 
try, genericity and regularity. In the classic example of 
the Ising model for ferromagnets, the expansion of the 
free energy as a function of the magnetisation m reads, 
in the homogeneous case and subject to external field h: 
F[m] = Fq — rah + |(T — T c )m 2 + Jjm 4 + from which 
a certain number of well known mean-field properties can 
be derived. Including the first gradient corrections in the 
inhomogencous case also allows one to show that close to 
T c , the divergence of the uniform susceptibility x(q = 0) is 
accompanied by the divergence of the correlation length £, 
over which magnetisation fluctuations are correlated. The 
Landau construction can falter in three distinct ways: 



1. Higher order terms, neglected in the expansion of 
F[m] as a series of m, could qualitatively change the 
above predictions (structural instability). This hap- 
pens, for example, close to a multicritical point where 
g(T c ) = 0. But if the transition remains second 
order, higher order terms are truly negligible when 
e = \T — T c \ — > and the predictions are universal. 

2. The non-linear feedback of spatial fluctuations on the 
divergence of the susceptibility can change all the crit- 
ical exponents when the dimension of space is smaller 
than d u = 4 in the case of the Ising model. For 
d > d u , on the other hand, one can prove that the 
low-q behaviour of x{l) is (close to the critical point) 
identical to that predicted by Landau's theory. 

3. Non-perturbative effects can wipe out the transition. 
This is the case for example of the spinodal transi- 
tion: the system is not able to reach the critical point 
because of nucleation, which is an activated process. 

Even the analogue of point (1) is difficult in the case of 
MCT; the basic reason being that the order parameter is 
a not a scalar, but it is a time dependent function C(k, t). 
The proof that MCT is structurally stable with respect to 
the addition of higher order terms is already quite complex 
and this will be the scope of the present paper. Once this 
is achieved, one should still worry about points (2) and (3) 
above. As already mentioned, it was recently realized that 
diverging fluctuations and an upper critical dimension d u 
also exist for MCT [16-18] (see also [19, 20] for earlier 
insights). In order to complete our proof that MCT is 
stable, one should prove that spatial fluctuations can be 
safely neglected in d > d u and understand how close one 
can get to the critical point before non-perturbative (ac- 
tivated) effects impair the transition. We will completely 
disregard these issues in the present paper, and focus only 
on point (1). 

The case of the glass transition is quite different from 
standard critical phenomena. Several physical and for- 
mal problems prevent a direct analogy. The glass transi- 
tion seems to be a purely dynamical phenomenon: simple 
static, thcrmodynamical properties do not present any pe- 
culiarities as the liquid freezes into a glass 1 . The above 
Landau construction simply does not make sense in the ab- 
sence of the clear analogy of the free energy. This means 
that the order parameter in glasses cannot be a one point 
function (such as the magnetisation) but, instead, it is 
likely to be a two point dynamic correlation. The slowing 
down of the dynamics in glasses is found to be accompa- 
nied by the appearance of a plateau value /k in the relax- 
ation pattern of the dynamical structure factor C (k, t) . 
Since the appearance of a plateau coincides with increas- 
ing time scales, one expects that within a very long time 

1 It has however been argued that highly non-trivial static cor- 
relations, called point to set correlations, increase approaching the 
glass transition [21]. 
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interval (to be specified), the correlation function can be 
approximately written as: 

C(k,r)«/ k 5 k + r5C(k,r), 5C(k,r)<l. (1) 

The idea underlying our construction of a Landau theory 
for glasses is to consider <5C(k, r) as the analogue of the or- 
der parameter m and construct a general, structurally sta- 
ble, dynamical equation for <5G(k, r). A way to construct 
such an equation is to start from the exact dynamical evo- 
lution for G and the response function R that can be de- 
rived in the framework of various dynamical field theories, 
for example based on Dean's equation [22] for Brownian 
dynamics, on Fluctuating Hydrodynamics for Newtonian 
dynamics [4] (see also [5-8]), or on Langevin equations 
for p-spin models [9]. Although these theories give very 
different sets of equations, they can all be reduced to the 
following single equation in the ergodic region: 

r 

d T C{k,r)+TC(k,T) + J duE(k,T-u)d u C(k,u) = (2) 
o 

with initial condition C(k, 0) = Sk- The self-energy 
E(k, r) (or memory kernel in MCT terminology) is given 
by the sum of 2-particle irreducible (2-PI) diagrams built 
with G and R lines, see e.g. [6]. We do not specify the 
details of the field theory underlying this equation, nor 
the Feynman rules for the diagrams contributing to E: 
we just need that such a theory exists. We also stay in 
the high-T region, so that the system is at equilibrium: 
both G and R are then time translation invariant and 
the Fluctuation-Dissipation theorem holds at a diagram- 
matic level: T R(k, r) = — d T C(k, r). In this case the 
self-energy is a functional of the correlation function only. 
Eq. (2) has exactly the structure of the standard MCT 
equation for liquids, although there is no well defined pre- 
scription to build a consistent approximation for E; in this 
sense MCT is rather arbitrary and difficult to improve 
upon in a systematic way. The standard MCT results 
correponds to a self-consistent 1-loop approximation for 
E, E(k, t - a) = J $^V(k, p)G(k - p,t - s)C(p, t - a), 
where V(k, p) is an effective vertex. 2 When generalized 
to higher order diagrams, an important difficulty emerges: 
the non-locality in time of the corrections. Our goal is 
to prove that the main results of the standard MCT (or 
1-loop) approximation still hold. The proof is done as 
for usual theories. We start from some conjectures about 
the critical properties, such as the nature of the order pa- 
rameter and its critical properties, that are motivated by 
experimental and numerical findings. Then we show that 
they result from a Landau-like expansion, which allows 
one to assess their universal character and to fix the value 
of the critical exponents. We therefore assume that the 
order parameter is the dynamical correlation function and 
that displays the following features: 

2 Note however that there are complications related to fluctuation 
dissipation relations [5-8] . 
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Fig. 1: A set of relaxation curves C(r) close to the MCT tran- 
sition for the so called schematic model of MCT [3] and several 
values of A, which quantifies the deviation from the transition. 
As A — > 1 a shoulder develops in the relaxation pattern corre- 
sponding to the emerging /3-relaxation regime. 

• There is structural arrest: below some temperature 
T d , lim C(k, t) = / k 5k with / k > 0. 

T — >00 

• When e — (T — T d )/T d <C 1, the correlation function 
exhibits a two-step pattern with three well separated 
characteristic time scales, see Fig 1. We assume that 
there exists a diverging time scale 773 (e) where the dif- 
ference <5C(k, t) between G and the plateau value is 
small, of order r(e) <C 1. More precisely, the corre- 
lation function decay is decomposed into: (i) a short 
time regime, r ~ tq, where C(k,r) = Co(k, r), with 
C (k,T > 1) -> / k Sk: (ii) a /3-regime, r = sr^(e) 
with a = 0(1): 5C(k, S T fi ) = r(e)S k (l - / k ) 2 G(k, a); 
(iii) an a-regime, t = s'r Q (e) with ,s' = O(l): 
C(k, s'r Q ) = G Q (k, s') describes the final fall off of 
the relaxation. 

We assume further (and justify later) that the func- 
tion G(k, s) can itself be expanded in powers of r(e): 
G(k,s) = £r=i r "~ lG «( k > s )- 3 All functions G n are 
a priori singular at s — and s = 00, reflecting the 
fact that the behaviour of C must match the short time 
regime and a regime, where the deviation from the plateau 
ceases to be small. A crucial remark for the follow- 
ing is that any function G n will appear with a prefac- 
tor r n (e). These hypotheses turn out to be sufficient to 
generalize the MCT results in the /3-regime. First, it is 
clear that the above expansion of G(k, s) generates a sim- 
ilar expansion of the self-energy E[C] in the /3-regime: 
E(k, s) = YlnLi 7 ' n ~ 1 '£'n{k, s), where the E„ do not de- 
pend on e, but are some functional of C(p,t). The most 
generic functional form for E a priori includes contribu- 
tions from all three regimes: 

E(k, s) = E[{C7 (p, aV/j), r(e)G(p, s'),C a (p, s'r /r a )}}, 
(3) 

3 One could have also some regular (in e) contributions. However, 
as we shall show the only two possible values of r(e) are ^/t or e. As a 
a consequence, regular contributions will be automatically contained 
in the expansion. 
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but since the £„ should not depend on e, general argu- 
ments can be used to restrict the actual functional form of 
S n . Note that we have used the notation s' to stress that 
this equation is a functional relation which is non-local in 
time. Also, even if we had assumed that SC(k, srg) only 
contains a single term of the order of r(e) then we would 
have generated corrections of all orders in r(e) anyway. 
The reason is that the self-energy will contain all orders 
in r(e) as it can be found by expanding the above equa- 
tion to all order in r(e)G(p, s); this will feed back, via the 
Schwinger-Dyson equations, on G(p, s) itself. 

We now illustrate how this works for the lowest order 
terms T, , £1 and £ 2 . As we shall see higher orders are in 
fact irrelevant for our purpose. Clearly, the zeroth order 
term S (k, s) can only be a function of the wavevector k 
since in the limit e — > time scales separate: rg — > oo and 
Tp/T a — ► and therefore the previous equation implies 
that in S (k, s) all dependence on s drops out. The first 
order contribution £i(k, s) must read: 



oo 

Ei(k,s) = J du J K 1 (k,p;s,u)G 1 (p, 



u) 



(4) 



where, henceforth, we shall use the notation J p = J ( 27r P 3 ■ 
Any other combination containing some G n gives an ex- 
tra factor r"(e) and thus it corresponds to a higher order 
contribution. In the original time variables, the kernel K\ 
must have some regular shape with a span fixed by the mi- 
croscopic time scale. Therefore, in the rescaled variables 
u, s, K\ must be local inu-s and, to lowest order, is a 
S function; higher derivatives of the S function correspond 
to corrections smaller by at least a factor 10/173 which, as 
we shall see, turn out to be negligible even at order r 2 . 
Therefore, S^k, s) = J j^K^k, p)Gi(p, s). 

The second order term has a richer structure. First, 
there is a term similar to the first order one with G 2 in- 
stead of d: J ^K 2 (k,p)G 2 (p,s). But since the kernel 
K 2 is obtained, as K\, from the first derivative of the self 
energy with respect to r(e)G(p, s), one finds K 2 — K\. 
Second, there are terms quadratic in G\. 



only the combination u + v — s enters comes from causality 
and the separation of time scales. The full justification of 
the above form and other technical details 4 are presented 
in [23]. The first local term (in s), Ku y £, is like the usual 
MCT contribution, but the other two terms do not appear 
within standard MCT. The third term actually reduces to 
the second one plus local terms via integration by parts. 

This expansion is the main result in the construction 
of the Landau theory. The zeroth order equation in r(e) 
fixes the non-crgodic parameter such that jzrfe = S (k), 
as in standard MCT. Substituting the expansion of £ up to 
second order in r into (2) and using the expansion of C in 
the (3 timescale, one finally obtains, for T = T<j(l + e) and 
in Laplace space within the /3-regime (we have dropped 
zeroth order, as discussed above) : 



T d {l + e)G 1 {\i,z)- f K 1 (k,p)G 1 (p,z) 
Jp 



+ 



T d (l + e)G 2 (k, z) 
Tdf]<e 



K 1 (k,p)G 2 (p,z) 



+ 



+ r z T d {l + e){\- f k )zGi(k,z) = 



+ 



+ 



z(l - /k) 

r 2 j K 2 {k,p)G 2 {p,z) 
Jp 



(5) 



f f /i:ii i< (k,ki,k 2 )r[Gi(ki,r)Gi(k 2 ,T)](2 ! ) 

Jki Jk 2 

: / / ifii,„£(k,ki,k 2 )zGi(ki,z)Gi(k 2 ,20, 

Jki Jk 2 



Identifying the coefficients order by order produces a se- 
ries of equations. The first order fixes the yet unknown 
function r(e): the expansion (5) only contains terms with 
integer powers of e. They should be matched with powers 
of r(e). Inspection of (5) shows that there are two possi- 
bilities 5 : either r — e, or r = \fi. The first choice yields a 
time independent solution for G\ which is in contradiction 
with our hypothesis of a two-step relaxation with diverg- 
ing time scales. Hence r(e) = ^/e, precisely as for usual 
MCT, or at 1-loop order. This follows from the presence 
of a regular in T term in (2). 

The equation to order r = ^[i now reads: 



jT ^ J du J dvK u {d,k 1 M;s,u,v)G 1 (k 1 ,u)G 1 (k 2 ,v) T d Gi(k, z) = j -^K^k, p)Gi(p, z) 



(6) 



For the same reasons outlined above, time dependence of 
ifn(q, ki, k 2 ; s, u, v) is composed of 5-functions and their 
derivatives. Some thinking about the underlying diagram- 
matic structure of the theory allows one to be convinced 
that the general structure of Kn is, to leading order: 

.£rii(s, u, v) = K\\^8{s — u)8{s — v)+ 
Ku, n eS( u + v - s )(<9« + d v) + K 11M (d u + d v )5(u + v - s) 

where to simplify the notation we have dropped all wave- 
vector dependence in the above equation. The fact that 



This is the standard eigenvalue problem found within 
MCT that fixes the value of the critical temperature T d - 
It constrains G\ to be a product of wave-vector depen- 
dent and time dependent amplitudes, thus reproducing 



4 In particular one could think that the separation of timescales 
not only leads to delta function terms but also to constants. The lat- 
ter are absent. This can be shown using the general field theoretical 
expression of the self energy, see [23] . 

5 Actually, there are other possibilities that correspond to higher 
order MCT singularities, which have been called A„ [27]. In a usual 
Landau theory these correspond to tricritical, or even higher order, 
critical points. We will neglect them here since they require some 
fine tuning of the coupling constants. 
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the well-know MCT "factorization property": Gi(k, z) = 
g(z)Hi(k), where Hi is the right eigenvector of K\ with 
largest eigenvalue A = T^. At this order however the 
scaling function g(z) remains unfixed. The second order 
equation is trickier: 



T d G 2 (k, «)- jT/f 1 (k,p)G 2 (p,«) = 



Tdfk 



-T d (l-f u )zg 2 (z)H 2 (k) + 



z(l - /k) 

+ f f ^i 1 ^(k,k 1 ,k 2 )£[ ff 2 ](^)F 1 (k 1 )F 1 (k 2 ) 
+zg 2 (z) ( ( ^ 11 ,„ < (k,ki,k 2 )Hi(ki)H 1 (k 2 ) 



(7) 



ki Jk 



Following [24], we now multiply (7) by -ffi(k) and integrate 
over k. The G 2 part of the equation vanishes and the 
remainder yields an equation on g{z). After some algebra 
and a proper rescaling of z and g one finds: 



-+ f {z )=C[g 2 ](z) 

Z A 



(8) 



where A is a constant that includes a non-local contri- 
bution as compared to MCT. But the structure of the 
equation on the scaling function g is exactly the same 
as in standard MCT. The properties of solution are well 
known: g has a singular power law asymptotics at z — > oo: 
g{z) ~ z a_1 and z — > 0: g(z) ~ The small time 

exponent a and long time exponent 6 characterize the de- 
cay of 5C(k,r) to and away from the plateau /k- The 
exponents a and b are related by the famous equation: 



r 2 (i- a ) _ r 2 (i + 6) 

r(l - 2a) ~~ r(l + 26) 



= A, 



(9) 



which is a genuinely non-trivial and clear-cut prediction 
of MCT that constrains the range of values of a and b 
to: < a < 1/2 and < b < 1, in good agreement with 
experimental and numerical results. We have thus found 
that this relation has a much broader degree of validity and 
survives the introduction of an arbitrary number of loop 
corrections. The values of a and b are however different 
from the standard MCT (or 1-loop) result, but as alluded 
to above, the parameter A is usually taken as an adjustable 
parameter anyway. 

The fact that the form of the scaling function g is the 
same as at 1-loop (MCT) has two consequences. First, it 
fixes the functional dependence of the time scales exactly 
as for MCT: 



-l/2a. 



(10) 



7 = l/2a + 1/26. This is clear from the matching of 
C(k, sTfj) at both ends of the (3 regime. Second, it can be 
used to show that some superficial divergencies encoun- 
tered in the calculation are in fact innocuous (see [23] for 
more details). Note that the only extra contribution that 



appear in the generic case to second order in g(s), namely 
the non-local term proportional to J dug' (u)g(s — u), does 
not modify the basic MCT equation, Eq. (8). 

The conclusion, which is the main result of this work, is 
that although T d and A are modified by taking into account 
corrections to MCT, the relation between the exponent a 
and 6, the square root singularity as well as the scaling 
function g are truly universal properties. This universality 
with respect to higher order local (in time) corrections was 
of course already shown by Gotze long ago; here we have 
proven that this result is robust with respect to general 
non-local corrections as well, and suggesting that MCT 
has the status of a Landau theory of the glass transition 6 . 

The above schematic arguments can be made precise 
within the context of specialized model. We have in par- 
ticular studied in full details the finite N corrections to 
mean-field 3-spin glass model, where the structure of the 
perturbation theory can be used to check that the above 
conclusions hold in that case, see [23] for details. 

It was recently understood how MCT equations should 
be generalized in the presence of spatial inhomogencities, 
where the correlation function C can be space dependent: 
C(k, r; r), where r is the average of the two points f\, r 2 
between which the correlation is computed, and k is the 
Fourier vector corresponding to f\ — r 2 . When wavelength 
of inhomogencities is large, one can establish a gradi- 
ent expansion of the MCT equations. In the schematic 
limit where all dependence on k is discarded, the self- 
energy reads, to the lowest order: E[C](s) = C(r, s) 2 + 
w\C(r, s)V 2 C(r, s) + ui 2 VC(r, s) ■ VC(r, s), where Wi and 
u>2 are some coefficients [17]. As mentioned in the intro- 
duction, these gradient terms are very important because 
they show how the MCT transition is in fact associated 
with a diverging correlation length, which corresponds to 
the scale over which a localized perturbation affects the 
surrounding dynamics [17]. The long-ranged critical fluc- 
tuations renormalizc the value of the MCT exponents in 
d < d u = 8 [18]. The above analysis, which was done 
in the homogeneous limit V — > 0, should be repeated in 
the inhomogencous case to complete our proof. We ex- 
pect that the same conclusion will hold, namely that the 
results about dynamical correlation obtained within inho- 
mogeneous MCT [17] are stable against the addition of 
higher order corrections. 

In conclusion, we have shown that MCT, which de- 
scribes a specific slowing down mechanism through the 
progressive disappearance of unstable directions, has the 
status of a Landau theory and is therefore expected 
to make generic predictions, albeit polluted by acti- 
vated events and critical fluctuations in finite dimensions. 
The interplay between critical fluctuations and activated 
events when d < d u , and the crossover to low temperature 
dynamics is still largely an exciting open problem [28]. 



6 As a consequence, it has a very different status compared to 
Mode-Coupling theories developed to compute critical exponents be- 
yond mean field theory for standard critical phenomena. 
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Note also that even for the exact MCT equations, the 
critical region where the asymptotic scaling predictions 
are valid is unusually narrow [25,26]. It would be in- 
teresting to generalize our Landau approach to the aging 
regime and show what are the truly universal properties 
of the mean-field and MCT-like description of the aging 
dynamics [9]. 

In constructing the Landau theory, we have assumed 
that the freezing transition is discontinuous, with a finite 
value /k of the plateau at the transition. A viable al- 
ternative is of course that of a continuous transition of 
the spin-glass type, which leads to a completely different 
phenomenology. This raises the question of the possible 
realization of this second scenario in the context of super- 
cooled liquids. All short-range interacting glasses seem 
to be characterized by rather small Lindemann parame- 
ters at the transition, meaning that it is hard to main- 
tain any kind of amorphous long range order when indi- 
vidual molecules move substantially, and that the glass 
transition is therefore discontinuous [29]. This argument 
suggests that continuous glasses can only exist for long- 
ranged interacting particles or quantum systems. In the 
quantum case, it is imaginable that amorphous density 
waves can indeed form with a vanishing modulation am- 
plitude (sec [30]). It would be very interesting to find 
experimental realizations of such a scenario. 
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